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The fixed point actions for Wilson and staggered lattice fermions are determined by iterating renormalization 
group transformations. In both cases a line of fixed points is found. Some points have very local fixed point 
actions. They can be used to construct perfect lattice actions for asymptotically free fermionic theories like QCD 
or the Gross-Neveu model. The local fixed point actions for Wilson fermions break chiral symmetry, while in the 
staggered case the remnant f/(l)e ® t/(l)o symmetry is preserved. In addition, for Wilson fermions a nonlocal 
fixed point is found that corresponds to free chiral fermions. The vicinity of this fixed point is studied in the 
Gross-Neveu model using perturbation theory. 
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1. Wilson Fermions 

The continuum limit of a lattice field theory 
is defined at a fixed point of the renormaliza- 
tion group. The lattice models on a renormal- 
ized trajectory emanating from the fixed point 
are free of cut-off effects and hence have perfect 
lattice actions. Recently, Hasenfratz and Nieder- 
mayer realized that perfect actions can be con- 
structed explicitly for asymptotically free theo- 
ries 1^] . In addition, in the 2-d nonlinear cr-model 
the renormalization group transformation can be 
optimized such that the fixed point action is ex- 
tremely local. This is essential for numerical sim- 
ulations. The question arises if fixed point actions 
for fermionic theories are local as well Hj. Since 
the fixed point of an asymptotically free theory is 
close to the Gaussian fixed point this question can 
be studied perturbatively, to lowest order even in 
the free theory. The corresponding calculation for 
a free scalar field was done long ago by Bell and 
Wilson §. _ 

Let us consider free Wilson fermion fields ^ and 
^ with the action S'[5', ^] on a hypercubic lattice 
A, which is then blocked to a lattice A' of doubled 
lattice spacing. Then each point x' £ A' corre- 
sponds to a hypercubic block of 2'^ points x E A 
and each point x belongs to exactly one block x' 
(we denote this by x G x'). The block transfor- 
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Figure 1. Blocking of a 2-d lattice. 



mation is illustrated in fig.l. On the blocked lat- 
tice we define new fermion fields 5" and with 
an effective action ,^'], and we perform a 
renormalization group step that leaves the parti- 
tion function unchanged 

exp(-5'[l'',1'']) = j P*X>*exp(-S'[l',*]) x 



2 



V^V^VrjVrj cxp(-S'[4', *]) 



exp{ [(4-1, - ^ 5^ *.)??x' + 



(1) 



In the last step auxiliary fields r}, rj have been in- 
troduced. For a = 00 the renormalization group 
transformation corresponds to a chirally invari- 
ant Grassmann 5-function, whereas for finite a 
the ri,j.,rix' term breaks chiral symmetry explic- 
itly. The parameter b is needed to renormalize 
the blocked fermion field. 

Starting from a point on the critical surface and 
iterating the renormalization group transforma- 
tion one generates actions on coarser and coarser 
lattices that may converge to a fixed point action 
S'*[*, *]. To find the fixed point we make the 
ansatz 

S[^, *] = « ^ - y)*x7M*y + 
x,y 

^X{x-y)^^^y. (2) 

x,y 

We go to momentum space and integrate out 5", 
and T], T] to obtain the blocked action S' , . 
At this point it is convenient to introduce 

""^^^^ " p(fc)2 + A(fc)2 ' ^^^^ = i l(fc)2 • 

After n renormalization group steps one finds 
a(r)(fc) = (6V2'^r$:a,(m^)x 



sin(fcj,/2) 



n 



/3(")(fc) = (6V2rE/?(^^)x (4) 



2"sin((A;^ -F27r/^)/2"+i)y ' 
2" 



n 



sin(fc^/2) 



1 - (672^)" 

^2"sin((A;^-h27r/^)/2"+i)y/ ^ a{l - 672'') 



Figure 2. The function ip\{z) for the fixed point 
action in coordinate space at a = 4. 



The summation is over integer component vectors 
I with e {1,2, 3,. ..,2"}. In the limit n ^ oo 
one approaches a nontrivial fixed point only if 



(5) 



The exponent (rf — l)/2 is the dimension of a free 

fermion field in d dimensions, and the blocked 
fermion field gets renormalized appropriately. At 
the fixed point one finds 



2/a. 



2ttI,. 



sin(fcy/2) 
{k + 2-Kiy J- J- \Kj2 + ^ 



n 



(6) 



Note that the result is independent of the param- 
eter r in the original action, as it should since the 
Wilson term is irrelevant. One finds a whole line 
of fixed points parametrized by a, and one may 
optimize a such that the fixed point action is as 
local as possible. In d = 1 the fixed point action 
has only nearest neighbor couplings when a = 4. 
This value is close to optimal also in d = 2. In 
fig. 2 the corresponding function p\{z) is displayed 
in coordinate space. Like in the scalar field case 
it is extremely local and very promising for nu- 
merical simulations of perfect fermion actions. 

For finite a both the renormalization group 
transformations and the fixed point actions break 
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Figure 3. The function pl{k) for the chirally in- 
variant nonlocal fixed point action. 



chiral symmetry and the fixed point actions are 
local. For a = oo, on the other hand, the renor- 
malization group transformation is chirally sym- 
metric. As a consequence, there exists a chirally 
invariant fixed point with jO*(fc) = a* (fc)/a*(fc)^, 
A*(fc) = 0, for which the Nielsen-Niniomiya theo- 
rem Q would suggest that the fermion spectrum 
is doubled. Fortunately, this is not the case be- 
cause the corresponding fixed point action is non- 
local and the theorem does not apply. Hence, the 
resulting continuum theory describes free chiral 
fermions. The function pl{k) is shown in fig. 3 for 
d = 2. The poles at the boundary of the Brillouin 
zone give rise to the action's nonlocality. SLAC 
fermions |^ also have a nonlocal action which is 
known to cause problems in perturbation theory 
|6). For fixed point fermions, however, the non- 
locality should be acceptable, because it arises 
naturally due the integration over the high mo- 
mentum modes of the fermion field. 

The real issue is to find a fixed point that de- 
scribes interacting chiral fermions. In particular, 
one should switch on a small gauge coupling and 



investigate the vicinity of the nonlocal fixed point 
of the free theory. Some time ago Rebbi Q sug- 
gested a lattice action for chiral fermions with a 
similar nonlocality, which unfortunately suffered 
from spurious ghost states ||,^. In fact, there is a 
no-go theorem due to Pelissetto [|l^ that excludes 
lattice chiral fermion constructions using certain 
nonlocal actions. Still, we are optimistic that in 
case of the nonlocal fixed point action the renor- 
malization group is clever enough to circumvent 
these arguments, although they certainly apply 
to generic man-made actions. 

So far we have studied the vicinity of the non- 
local fixed point in the Gross-Neveu model. The 
4-Fermi interaction is linearized by an auxiliary 
scalar field $ with a Yukawa coupling y. In 
momentum space the Yukawa coupling takes the 
form 



d2fcid2fc2^(fci)a(fci, A:2)^(fc2)$(-fci - fc2). (7) 



The coupling a{ki, is a matrix in Dirac space. 
Now the renormalization group is iterated by also 
blocking the field <& to coarser and coarser lat- 
tices. At the fixed point of the coupled system 
one finds to leading order in y 

a*{ki,k2) = y p^ikih^ ^ fci,. + 27r?i^ 
k2p + 2nl2p 



(fci + 27r/i; 



Til- X 



(fc2 

n 



27r?2)2 

sin((fciA 



7p P*aik2ha X 

k2x)/2) 



(fciA -I- 2t:Ii) 
sm{kix/2) 



kax + 2^l2x)l2 
sin(fc2A/2) 



{kix + 2^hx)/2 (fc2A + 2^?2A)/2' 



(8) 



Before attacking the gauge theory we will use this 
vertex function to check if ghost states spoil the 
continuum limit in this nonlocal chirally invari- 
ant formulation of the Gross-Neveu model. This 
investigation is in progress. 

2. Staggered Fermions 

An alternative formulation of lattice fermions 
uses staggered fermion fields. This has the advan- 
tage that the cut-off theory has a remnant chiral 
symmetry. The class of lattice fermions which 



4 



is described by the ansatz (g) does not include 
staggered fermions, since their action is invariant 
only under translations by an even number of lat- 
tice spacings. Staggered fermions can be treated 
analogously, but the renormalization group trans- 
formations are different and new properties of the 
fixed point actions arise. 

As usual we view the staggered fcrmion vari- 
ables Xx, Xx fields with pseudoflavors i = I + 
ni + 2n2 + ... + 2'^~^nd {n^ G {0, 1}) located at the 
2'' corners of the hypercubes centered at points x 
which form a lattice with spacing 2. Staggered 
fermions have various symmetries, among them 
a U{l)e ® U{\)o remainder of chiral invariance, 
an analog of charge conjugation Co, and shift 
symmetries containing translations by one lattice 
spacing. To be consistent with these symmetries 
the renormalization group transformation has to 
be modified. At the end, one wants to reconstruct 
Dirac spinors from the pseudoflavors. Therefore 
it is important not to mix the corners of the hy- 
percubes in the renormalization group transfor- 
mation. Kalkreuter, Mack and Speh have pro- 
posed a suitable blocking scheme with block fac- 
tor 3 (in general it must be odd) . One among 
S'' block centers x remains a block center after the 
renormalization group step. We denote it by x' . 
The x' form a lattice of spacing 6. Each pseud- 
oflavor builds its individual block variable, and 
every on the fine lattice contributes to exactly 
one x^'/ on the coarse lattice. The block transfor- 
mation is illustrated in fig. 4. First we apply a 
5-function renormalization group transformation 
using this blocking scheme 

exp(-5'[x',x']) = j VxDxeM-S[x,x]) x 

\{^^^' - |i E ^)^^^' - |j E 

x' ,i x^x' xGx' 

= J VxVxVfjV7jeM-S[x,x]) x 

«-p{E E[(x^'-|iEx;K' + 

x' ^A' i x£x' 

nx'ix2'-^T.^^}- (9) 

xGx' 

For each pseudoflavor we have defined coarse 
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Figure 4. Blocking of a 2-d lattice consistent with 
the staggered symmetries. 



grained auxiliary Grassmann fields ry* / , ry^/ . 
Clearly this is analogous to the chirally invariant 
transformation eq.(^) for Wilson fermions with 
a = oo. 

Now we look for a general ansatz, analogous to 
(^), which describes the action after a number of 
renormalization group steps 

S[x,x]^i^^xlPtj{x-y)xl- (10) 

Here p depends only on the distance x—y between 
two hypercube centers, because the action is in- 
variant against translations by an even number of 
lattice spacings. 

The standard action as well as the renormal- 
ization group transformation (^ respect various 
symmetries, which therefore also hold for (p^). 
This specifies the form of the matrix p as fol- 
lows: the U{l)e®U{l)o symmetry implies pij ^ 
only if the pseudoflavors i and j belong to dif- 
ferent sublattices (for one pseudoflavor the sum 
of lattice point coordinates is even and for the 
other one it is odd), and Co invariance requires 
Pij{x — y) = —pji{y — x). Shifts by one lattice 
spacing yield further relations among the remain- 
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ing degrees of freedom; they relate sets of 2'^~^ 
matrix elements. As for Wilson fermions it is nat- 
ural to go to momentum space and to work with 
the inverse matrix a(fc) = p{k)~^. In particular, 
for d = 2 the various symmetries imply 



a{k) 





ai{k) 
a2{k) 




ai{k) 



-a2{k) 



ceiik) 



ai{k) 





-a2{k) 
ai{k) 




(11) 



In higher dimensions also couplings along cer- 
tain space diagonals arc permitted by the sym- 
metries; however, they are absent in the stan- 
dard action and it turns out that the renormaliza- 
tion group transformation docs not activate them 
either. Hence, in general one can parametrize 
the inverse of p{k) by d functions a^(A;) with 



a, 



After n renormalization group steps we obtain 



Figure 5. The function ipl2{z) for the staggered 
fermion fixed point in coordinate space (d = 2). 



n 



sin(fci//2) 



3" sin((fc^-|-27rZ^)/3"2) 



(12) 



with G {1,2,3, 3"}. The sign factor (-l)'*- 
arises because certain terms are antiperiodic with 
respect to fc^ . As for Wilson fermions a nontrivial 
fixed point is reached only for one particular value 
of b; in this case b = 3^"^"^'/^ in accordance with 
dimensional considerations. 

It turns out that the U{l)e ® U{l)o invariant 
fixed point action is local as one infers from figs. 5 
and 6. This is an important qualitative differ- 
ence to the chirally invariant nonlocal fixed point 
for Wilson fermions. In the staggered case there 
is no contradiction to the Nielsen-Ninomiya the- 
orem since the U (l)e d) U (l)o symmetry does not 
imply full chiral invariance of the Dirac spinors 
that one can reconstruct. Adding a mass term 
of the auxiliary fields (as it was done to optimize 
the locality of the fixed point action for Wilson 
fermions) would break the remnant chiral sym- 
metry and would therefore destroy an essential 
advantage of the staggered fermion formulation. 
Still, we can improve the locality without loss of 



Figure 6. The function Pi2(^) the staggered 
fermion fixed point action (d = 2). 
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Figure 7. The function ipl2{z) for the optimally 
local fixed point {d = 2, a = 9/4). 



symmetry if we add a kinetic term of the auxil- 
iary fields, again suppressed by a parameter 1/a. 
At the fixed point one obtains 



lez'' 



n 



+ nliy 



■ sin(fc^/2). 



(13) 



As long as no gauge fields are present one can still 
show that the partition function remains invari- 
ant under renormalization group transformations. 
In d = 1 the fixed point action is an optimally lo- 
cal nearest neighbor interaction when a = 9/4. 
Numerically it turns out that this choice is also 
successful for d = 2 as one sees in fig. 7. However, 
if we switch on a gauge field this transformation is 
not allowed any more and we are restricted to the 
^-function transformation (a — oo) from above. 
In summary we repeat that the fixed point ac- 
tions we obtain for staggered fermions are local 
and U{l)e ® U{l)o invariant, and are hence well 
suited for numerical simulations. 



Next we intend to search for the fixed point 
action in the Gross-Neveu model using staggered 
fermions. Here the kinetic term of the auxiliary 
field can be used and the location of the critical 
surface is obvious. Since the Gross-Neveu model 
is asymptotically free, there is one (weakly) rel- 
evant direction (to lowest order it is marginal). 
When one determines this direction one obtains 
the analog of the line of classical perfect actions 
in the bosonic case Q. Then one may adopt the 
concept of Hasenfratz and Niedermayer and fol- 
low this (straight) line far enough to small correla- 
tion lengths where numerical simulations become 
feasible. If the simulated action is close to a per- 
fect quantum action one can extract continuum 
physics to a good approximation with only little 
numerical effort. However, in the fermionic case 
it remains to be seen if an (almost) perfect action 
can be simulated with the Hybrid-Monte-Carlo 
algorithm, or if one is restricted to less efficient 
numerical methods. 

We are indebted to P. Hasenfratz and F. Nie- 
dermayer for very interesting discussions. We also 
like to thank T. Kalkreuter for providing us with 
a copy of fig. 4, and A. Pelissetto for explaining 
his no-go theorem. 
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